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Abstract 


Since the introduction of the finite element approach, as a numerical solu- 
tion scheme for structural and solid mechanics applications, various for- 
mulation methodologies have been proposed. ‘These ways offer differ- 
ent advantages and shortcomings. Among these techniques, the stan- 
dard displacement-based approach has attracted more interest due to its 
straightforward scheme and generality. Investigators have proved that the 
other strategies, such as the force-based, hybrid, assumed stress, and as- 
sumed strain provides special advantages in comparison with the classic 
finite elements. For instance, the mentioned techniques are able to solve 
difficulties, like shear locking, shear parasitic error, mesh sensitivity, poor 
convergence, and rotational dependency. The main goal of this two-part 
study is to present a brief yet clear portrait of the basics and advantages of 
the direct strain-based method for development of high-performance plane 
finite elements. In this article, which is the first part of this study, assump- 
tions and the basics of this method are introduced. Then, a detailed review 
of all the existing strain-based membrane elements is presented. Although 
the strain formulation is applicable for different types of structures, most 
of the existing elements pertain to the plane structures. The second part 
of this study deals with the application and performance of the reviewed 
elements in the analysis of plane stress/strain problems. 


AMS subject classifications (2020): 74K15, 74G15. 


*Corresponding author 

Received 25 October 2020; revised 8 June 2021; accepted 9 June 2021 

Mohammad Rezaiee-Pajand 

Professor of Civil Engineering, School of Engineering, Ferdowsi University of Mashhad, 
Iran. e-mail: Rezaiee@Qum.ac.ir, Tel/fax: +98-51-38412912 


Nima Gharaei-Moghaddam 
PhD of Structural Engineering, School of Engineering, Ferdowsi University of Mashhad, 
Iran. e-mail: Nima.Gharaei@gmail.com, Tel: +98-915-1589342 


Mohammadreza Ramezani 

PhD Student of Structural Engineering, School of Engineering, Ferdowsi University 
of Mashhad, Iran. e-mail: Mohammadrezaramezanil994@gmail.com, Tel: +98-915- 
1076010 


437 


438 Rezaiee-Pajand, Gharaei-Moghaddam and Ramezani 


Keywords: Strain-based formulation; Higher-order strain field; Equilibrium 
condition; Numerical evaluation; Drilling degrees of freedom. 


1 Introduction 


Numerical methods are proved to be powerful and effective computational 
tools for analysis of complicated and practical engineering problems. Ac- 
cording to different features of scientific activities and in demand of spe- 
cial requirements, many diverse numerical techniques are developed in the 
past decades, such as the finite element method, finite difference technique, 
boundary element method, and discrete element approach. Each of these 
various numerical methods have their own advantages and shortcomings, but 
the finite element methods gain more popularity due to their strong math- 
ematical bases and inherent capabilities, which result in increasing applica- 
tion of this scheme in different fields of science and especially engineering 
fields [5, 6, 28, 39, 40, 42, 43, 78, 79]. Therefore, various formulation tech- 
niques are developed in the past decades and there are thousands of finite 
elements available for analysis of dissimilar types of problems and structures 
[8, 13, 21, 22, 41, 53]. 

Among the available approaches for finite element formulation, the most 
well-known and widely applicable one is the displacement-based technique. 
This method, which sometimes is called with different terms, such as, the 
classical or stiffness approach, is the first scheme that was used for the de- 
velopment of finite elements; see [80]. Clear and straightforward process and 
applicability to different types of problems and structures are the promi- 
nent advantages of the displacement-based formulation for structural and 
mechanical applications. However, this process has various shortcomings. 
For instance, inaccuracy and discontinuity of stresses, which are secondary 
parameters in stiffness approach, are a vital deficiency in structural applica- 
tions, where stress is a decisive parameter in the design practice; see [83, 52]. 
It should be noted that this problem can partly be solved by using higher 
order formulations, which leads to the application of internal nodes, and 
increases computational costs, and reduces the numerical efficiency of the 
analysis [23, 36, 54, 55]. Another common problem of displacement-based 
finite elements in various locking phenomena, such as, shear and membrane 
locking, which necessitate special treating, which sometimes requires con- 
siderable time and effort and reduces the efficiency of the method [54, 55]. 
Moreover, in severely nonlinear problems, the displacement-based elements 
usually necessitate utilization of very fine meshes, which is inappropriate 
from the efficiency standpoint [46, 47, 74]. To remedy the mentioned and 
other shortcomings of the displacement approach, other finite element for- 
mulations, such as the force-based, hybrid or mixed, assumed stress, and 
assumed strain has been developed. Fortunately, these new procedures have 
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their own advantages and shortcomings. For example, the force-based for- 
mulation performs very well in the linear and nonlinear analysis of frame 
structures and also provides an appropriate platform for the development of 
advanced frame elements [3, 34, 56, 57, 73, 76]. It is reminded that the force 
formulation approach is limited to skeletal structures and that its application 
for continuous structures is very difficult if not impossible. 

Although each method has different merits and limitations, some of the 
approaches have received more attention from researchers, while the others 
have remained less treated. One of these techniques that has received less 
attention despite its promising performance, is the strain-based or assumed 
strain approach [17, 25, 35, 38]. Therefore, the main purpose of this study is 
to introduce basics of the strain-based formulation. The strain formulation 
method can be classified in three distinct groups, namely free formulation, 
assumed natural strain method, and direct approach. 

The free formulation method is based on the kinematic decomposition, 
in which the element displacements are decomposed in two basic and higher- 
order parts [25, 38]. The basic part represents the rigid body motions and 
constant strain state, which is necessary for the convergence criteria. On the 
other hand, the higher-order terms improve accuracy of the finite element 
by establishing proper rank of the stiffness matrix. The conditions of force 
and energy orthogonality result in the algebraic formulation of the stiffness 
matrix, which satisfies the individual element test. It is noteworthy that 
the strain functions in this method are dependent on displacements one; see 
[19]. The second method for the development of strain-based elements is the 
assumed natural deviatoric strain (ANDES) approach [27]. In this technique, 
the independent strain function is applied, which includes basic and higher- 
order parts. The deviation of strain from the constant strain is represented 
by the higher-order part of the strain function. In this way, the higher-order 
part is selected so that the integral of higher-order strain through the element 
equals zero. Therefore, ANDES technique satisfies the individual element test 
[27]. 

The third method is the direct formulation approach, in which Taylor 
series for the strain field is used to approximate the strain field [63]. The 
resulting elements from this approach lead free from the shear locking and 
parasitic shear error. The strain states in this approach can be divided into 
rigid body motions, constant strain, and higher-order strain states. The rigid 
body modes and constant strain states guarantee the convergence of result- 
ing element. The higher-order terms have a parametric form that can be 
optimized to obtain efficient elements. The optimization is performed by en- 
forcing different optimal conditions. Some of these optimal conditions will 
be discussed in the coming sections. It is important to note that the opti- 
mization process of the finite element templates is a relatively complicated 
task that requires innovation [63]. 

As mentioned, the strain-based formulation itself can be categorized in 
three different subdivisions. In this study, only the direct method is covered. 
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After the presentation of the formulation basics, a detailed review of the 
existing assumed-strain membrane elements is presented. This article is the 
first part of a two-part study. In the second part, numerical comparison of 
the reviewed elements is performed. This article is organized in the following 
order: Section 2 presents basis steps for the development of the plane elements 
with the assumed strain approach. Other important optimal criteria for 
assume strain elements are introduced in Section 3. The existing triangular 
membrane elements are reviewed in Section 4, and Section 5 discusses the 
available quadrilateral plane elements. In Section 6, some of the accessible 
strain-based elements for the other types of structures are briefly introduced. 
Finally, Section 7 presents the concluding remarks of the article. 


2 Basics of the formulation 


In the assumed strain formulation, the element strain field is approximated 
by an assumed mathematical function. As mentioned, like any other form 
of the finite element formulation, there are different types of assumed strain 
formulation [4, 20, 26, 37]. However, in this study, only the direct strain- 
based formulation will be discussed [63]. The following formulation steps 
pertain to development of plane finite elements. Needless to say, the pre- 
sented process can be simply applied to the other types of elements with only 
slight modifications. The main reasons behind narrowing the scope of the 
present review are the availability of the relatively large number of existing 
strain-based elements developed based on the various types of strain-based 
formulation. These researches made the current study very lengthy. More- 
over, the growing interest from the structural analysis community toward 
the direct method in recent years, and also the valuable experiences of the 
authors in this field, which should be exposed to the younger analysts. 

In the case of plane problems, the strain field consists of three components, 
namely €z, €y, and Yzy. Based on the Taylor expansion, each function can 
be approximated by a polynomial function of arbitrary order. Therefore, the 
strain components are approximated as follows: 


2 
ex (2,9) = (fx)o + (Sx,2) 9% + (fx,y)o¥ + (Cx,22)q (3) + (Cx,2y)o (#u) + (Cx.vy)o (4) ++ 
2 2 
ey (20) = (eyo + Cue)o® + (Quulow + eu.e2)o (3) + ev.eu)o (ev) + (evuy)o (Ye) te. 
2 
yey (@, 9) = (Yey)o + (Yey,2) 5% + (Yey,y)o¥ + (Yay,er) 6 +) + (vay, zy), (ey) + (Yay, yy)o (+) Sosa 


(1) 
Here “,” indicates differentiating with respect to its following variable. More- 
over, the subscript “o” indicates the value of the strain gradient at the origin 
of the coordinate system. The coefficients with the subscript “o” are called 
strain states. Selection of diverse polynomials with a different number of 
terms and various orders results in finite elements with a different number 
of degrees of freedom, as well as, the specific properties. Due to the impor- 
tance of the constant strain states for the convergence of the resulting finite 
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element, the selection of the constant and linear terms for strain components 
is necessary. Selection of the higher-order terms would increase the accuracy 
and the convergence rate of the resulting finite element, but instead reduces 
its numerical efficiency. As it was mentioned previously, it is possible to ap- 
ply different optimal criteria, such as pure plain bending test, for improving 
performance of the resulting finite element. However, such optimal condi- 
tions are optional, and the only necessity is to include the constant terms 
in the opted assumed strain field. However, similar to the classical formu- 
lation, it is suggested not to give any priority to each of the coordinates, (x 
or y). In addition, if the analytical estimation of the strain field is available, 
then the terms of the approximated field can be selected based on the known 
analytical solution. 

Nevertheless, after choosing the desired terms, the assumed strain field 
can be optimized by applying any desired optimal condition [38, 63, 71]. 
The most common criteria are compatibility and equilibrium conditions [63, 
71]. The compatibility of the strain field is achieved provided that the next 
relationship exists between the strain components [71]: 


Cee Ob, _ Ovary 
Oy” Ox? Oxdy” 


(2) 


It is obvious that in the general case, for the satisfaction of the compatibility 
or any other criteria; it might be needed that some strain states be dependent 
on each other. Therefore, imposing the optimized condition results in the 
dependency of some strain states to each other and therefore, it reduces 
the number of independent strain states [27]. As it will be shown later, 
the number of independent strain states is equal to the number of required 
degrees of freedom for the element. 


The other conventional optimal condition is the equilibrium. It was 
proved that if the equilibrium equation is satisfied within an element, then it 
can be included in the Trefftz formulation [27]. The equation of equilibrium 
for the plane problems is defined as follows: 


Oox OTxy —_ 

Ox + Oy |r 0, (3) 
Ooy + Oey = 

Oy On Try 


where F, and Fy are the body forces in the z and y directions, respectively. 
Also, oz, Jy and Tzy are normal and shearing stresses, respectively. To 
rewrite the equilibrium equation in terms of strain, it is necessary to relate 


the stresses to the strains. For the plane problems, the coming relations 


connect stresses and strains: 


Oy = 2GEy + 


Oy = 2Ge, 4 


Try = Gray: 


Alég +€y)s 
A(Ex + Ey), 


(4) 
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In the previous equations, G and v are the shear modulus and Poisson’s ratio, 


respectively. Indeed A is called the Lame constant and is equal to Guo 


for the plane stress condition. In the case of plane strain, this constant is 


equal to Oe Here, F stands for modulus of elasticity. Substituting 


equations (4) in the equilibrium equation results in the following relations: 


a age (5) 
AGE + (8G + A) Ze + GH + Fy, =0. 


ies +) Sen 4 2x 4 GOtev + Fe, = 0, 


It is noteworthy that imposing the optimized conditions, such as; compati- 
bility and equilibrium, is not necessary steps for developing a plane element 
based on assumed strain approach, and as it will be shown in the coming sec- 
tion, there are finite elements, which do not consider these criteria [71, 72]. 
However, enforcing these conditions to the assumed strain field improves the 
performance of the resulting element. As mentioned previously, the inclusion 
of the optimized condition makes some of the strain states to be depended 
on the other ones. When the dependent strain states are determined, the 
assumed strain field is rewritten in terms of the independent ones. The next 
step is to calculate the associated displacement field. For this purpose, the 
strain-displacement formulas are utilized: 


Eg = gu 
Ey = a (6) 
Yay = 3 (5H + gv) 


In these relations, u and v are displacements in the xz and y directions, 
respectively. Based on these qualities, the displacements in x and y directions 
are derived by integrating normal strain components with respect to their 
associated coordinates: 

u(x,y) = f cada + fily), (7) 
v(e,y) = feydy + fale): 


Here, f; and fe are the results of integrating shear strain with respect to 
the coordinates and imposing the rigid body modes condition. In the case 
of plane problems, three rigid body modes exist in the displacement field, 
namely ug, Vs, and ro, signifying the rigid body displacements in x and y 
directions and the rigid body rotation, respectively. As it was mentioned, the 
existence of these terms is the necessary convergence condition. Accordingly, 
these modes are also counted among the independent strain states that can 
be arranged in a vector arrangement indicated by S. This vector is called 
strain state vector, which consists of the independent strain states. They 
are the coefficients of the strain field approximations in (1), as well as, the 
rigid body modes. The strain state vector is somehow equivalent to the 
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nodal displacement vector in the traditional displacement-based approach. 
By using the matrix notation, which is traditionally used in finite element 
formulation in structural engineering applications, it is possible to relate the 
displacement and strain fields to the strain state vector in the subsequent 
forms: 


U=N,S+U, (8) 
e=B,S+E. (9) 


In these equations, N, and B, represent the displacement and strain inter- 
polation matrices, respectively. The particular response of the displacement 
and strain fields, that is, U and €, depend on the body forces. The next 
relation can be established between the vectors of nodal displacements and 
the strain states: 

D=AS+D=D+D. (10) 


Here, D and D are the nodal displacement vectors and the displacements 
due to body forces. Also, A is a geometric matrix including of the nodal dis- 
placement interpolation matrices. It is possible to construct the subsequent 
relations between the displacement and strains’ fields of the element with the 
nodal displacement vector using (10) as follows: 


U=N,S+U=N,(A1D)+U=(NA")D+U=ND+U, (11) 
e=B,S+€=B,(A'D)+@=(B,A')D+F= BD +E. (12) 


Assuming the body forces to be negligible in comparison with the applied 
external loads, the strains and displacements due to body forces, € and U, 
are neglected. 


The final step of the formulation is to derive the element stiffness matrix 
and the nodal force vector. Among different approaches for this purpose, 
the minimization of the total potential energy is the most common approach. 
The total potential energy functional can be established as follows: 


Il = 5 fete Wa / UT F dV — D? Pont. (13) 


In this equation, P.., and F’ stand for the external nodal and the body 
forces, respectively. The element stiffens matrix and nodal force vector are 
derived by establishing the stationary of the following functional: 


24 (/ 7 Dm Bad) A-'D-AT ¢; Nv." Fav) Peat = KD—P=0. (14) 


oD 


Therefore, the element stiffness matrix and the nodal force vector are 
derived in the following form: 
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K=AT (/ B. DmB,dv) At1=ATK A“, (15) 


P= Poy, +AT (/ N," Fav) ; (16) 


where, D,,, is the material matrix: 


lv O 
Din = 5 yl 0}. (17) 
—~H' 100 t= 


It is noteworthy that for plane problems, which is the main subject of this 
study, the volume integral in the above-mentioned equations simply trans- 
forms into an area integral considering the constant unit thickness for the 
plane structures. 


3 Required optimal criteria 


As mentioned previously, in order to achieve an optimal and efficient for- 
mulation, it is possible to impose different criteria on the assumed strain 
field. Two of these criteria, compatibility and equilibrium, are defined in the 
previous section. In this section, two other required criteria for the optimal 
performance of assumed strain elements are defined. 


3.1 Pure bending test 


To achieve optimal performance in flexural behavior, Felippa [18, 26] utilized 
the pure bending test. In this experiment, an Euler—Bernoulli beam is dis- 
cretized by rectangular (or triangular) elements and loaded by the constant 
bending moments. To study in-plane bending along x and y axes, the two 
beams depicted in Figures 1 and 2 are considered. 


y y 
Mx 


Pipi 


Figure 1: Pure bending test in the zy plane 


The stored elastic energy due to the constant moments, M, and M,, in 
part of the beams meshed with one rectangular element (or two triangular 
elements) is derived according to the following relations: 


Review of the strain-based formulation for analysis of ..., Part I 445 


6aM2 
exact = x 1 
Usenet = (18) 
6bM? 
exact __ ¥y : 1 
Op Eath 2) 


i COO 


Pipe I 


beam cross section 


x 
EE 
<——> 

b 
Figure 2: Pure bending test in the zy plane 


Based on this test, an element can present bending behavior exactly, 
provided that it can compute precise elastic stored energy. In other words, 
the ratio of the energy calculated according to the results of element analysis 
to the exact stored elastic energy should be equal to 1. For this purpose, the 
element stored energy, based on finite element responses, is computed using 
the coming equations: 


1 
pone = 5 Phir K Dow: (20) 
1 
poe —_ 5 Poy K Dey: (21) 


In these relations, Dy, and Dp, are the nodal displacement vector of the 
element and K is the element stiffness matrix. Therefore, the previously 
mentioned flexural energy ratios are derived as follows: 


element 
—— U; 


"es = Ueract ’ (22) 
Uelement 
¥ 
Ty = TrexaGk (23) 
US t 


If rz or ry is equal to 1, the element passes the pure bending test and is 
able to represent exact bending behavior. If these ratios are greater or less 
than 1, then the element is over-stiff or over-flexible, respectively. Moreover, 
if the energy ratios are equal to 1 for any aspect ratio (a/b), then the element 
is optimal in bending behavior. Finally, the element suffers from shear locking 
on the condition that a >> b resultsinr, >> 1orifb >> aleadstor, >> 1. 

To study this test for the strain-based formulation, it is required to in- 
vestigate the exact strain fields of the elements under pure bending. In the 
following, only the case of M, will be discussed. The same reasoning goes 
for M,, as well. The exact stress components of the beam subjected to the 
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pure bending M, are as follows: 


Ox = — BA? 
oy =0, (24) 
Ty = 0. 


Therefore, the element stress field can be demonstrated by the linear function 
of y for oz in the subsequent form: 


Or =A, + Q2Y, 
Jy = 0, (25) 


Try = 0. 


Based on the Hook’s law, the corresponding strain field is as follows: 


Ee = G+ Gy =F. + Boy, 
Oy = et — ey = Bs + Bay, (26) 
Tay = 0. 


From the previous relations, it can be concluded that a strain-based element 
would surely pass the pure bending test, provided that the constant and 
linear terms for the normal strains are included in the assumed strain field. 


3.2 Rotational invariance 


Because the finite elements may be rotated with the different angles and 
be placed in various locations in the structural meshes used for the analysis 
of the diverse problems, their characteristics must not change due to the 
rotation. Such an element is called a rotational invariant. Assuming that 
the coordinate system zy is rotated to a new form, which is indicated by 
xz’y’. The displacements’ components in this new coordinate system can be 
related to those of the initial coordinate system by using the following simple 
transformations: 


u’ = ucos@ + vsin 8, (27) 
vu’ = —usin @ + vcos 6. (28) 


In this relation, @ is the rotation angle from zy system to the new x’y’ coor- 
dinates. Based on this relation, having the rigid body motions, ug and vp in 
the original coordinate, it is possible to calculate the correct value of up and 
up, in the x’y’ coordinates, by using equations (27) and (28). Regarding the 
rotational invariance property, it is required to consider the strain field with 
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a complete order. In other words, some incomplete interpolation polynomi- 
als produce strain states that are not invariant with rotation. Therefore, the 
rotational invariance can be guaranteed with the inclusion of all strain terms 
with a given order. For instance, the rotational mapping of a constant strain 
state is as follows: 


€,, = £0080 +€,sin?9 + 7xysin 6 cos 6 , (29) 
Ey = €,8in?9 +e,cos?0 —Yxysin @ cos @ , (30) 
Vey = (Ex — Ey) Sin 20 + yeycos 26 . (31) 


Based on these relations, a strain-based element can represent the constant 
strains with respect to any system of the coordinates, only on the condition 
that its formulation takes into account all three cases of the constant strain 
states. Although, the completeness of the assumed strain field guarantees 
rotational invariance of the element, the elements with incomplete strain 
fields are not necessarily rotational dependent. 


4 Triangular membrane elements 


To the author’s best knowledge, there are thirteen strain-based triangular el- 
ements formulated by using direct strain-based formulation. In this section, 
thirteen triangular membrane elements proposed based on the assumed-strain 
approach are briefly reviewed. These studies are arranged in historical order. 
Most of the available triangular element has similar geometry, a three-node 
triangle with three degrees of freedom at each node. However, in this con- 
figuration, the incorporation of the drilling degrees of freedom improves the 
performance of the element in bending analysis, but more recent works uti- 
lized different distribution of degrees of freedom provide better results. In 
addition, as it will be demonstrated, in the most of the existing elements, 
the equilibrium conditions are not imposed on the assumed strain fields. Ac- 
cording to the outcomes of the more recent elements, considering the equilib- 
rium conditions improve the accuracy and convergence rate of the suggested 
strain-based elements. More details about these elements are provided in the 
following descriptions. To facilitate understanding and reproducing of the 
elements by readers, the main details are presented in a table format. 


4.1 Sabir (1985) 


The first researcher that utilized the assumed strain approach to develop 
more powerful membrane elements is Sabir. In one of his early works, he [71] 
proposed a three-node triangular element with three-degrees of freedom at 
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each node. The assumed strain components satisfied the compatibility equa- 
tion, but the equilibrium equations were not satisfied. The drilling degree of 
freedom for the element was defined by the subsequent relation: 


1 fou Ou 


Details of the element formulation are presented in Table 1. 


4.2 Sabir and Sfendji (1995) 


Sabir and Sfendji [72] suggested a four-node triangular element by assumption 
of the linear normal strains and constant shear strain. The selected strain 
field satisfied the compatibility condition, but the equilibrium equations were 
not imposed on this strain field. Therefore, the strain state vector consisted 
of eight independent unknown coefficients. Three nodes were located at the 
vertexes of the triangle, and the fourth node was placed in the middle of one 
side. Each node had two translational degrees of freedom. Consequently, 
the element possessed eight degrees of freedom totally. This geometry made 
the element appropriate to be used as transitional element in finite element 
meshes. They compared the performance of their suggested element with 
standard displacement-based element by using few simple numerical problems 
[72]. Their attained results showed better performance of this triangular 
strain-based element. Table 2 presents details of this element. 


4.3 Tayeh (2003) 


In 2003, Tayeh [75] developed new strain-based elements for analysis of the 
plane structures. In contrast to the previous works by Sabir, he utilized 
higher-order terms and assumed an incomplete second-order field for the 
element (see Table 3). It is evident that some coefficients in this assumed 
strain field are common between different strain components. Tayeh provided 
no clear reason for this selection, but he stated indirectly that this strain field 
was selected in order to satisfy the compatibility condition. Similar to the 
element proposed by Sabir and Sfendji [72], the assumed strain field did not 
satisfy equilibrium equations. 


Review of the strain-based formulation for analysis of ..., Part I 449 


Table 1: Details of the triangular element proposed by Sabir [71] 


Ref | Properties Geometry 
Geometry 
Ex (@,y) = (Ex), ar (Ex,y) oY Tr (Ey,0) 9% 
Strain field Sy (2,4) = (Ey)g + eye) 8 + (Cea) 
= Yey (,Y) = ey) + Yay.) ot + (Yay,0) oY 
=. [Strain state vec-[S={us vo To (Ex), (Ey), Cey)o 
a tor (Ex) o (Ey,2)o (Veiga a 
Ww Nodal displace- D= {D3;i_2 D351 D3; D3;—2 
ment vector D35-1 D3; D3p_9 D3p_4 D3y}F 
000100ya O 
Strain interpola-|B,= |0O00010ya _ 0 
tion matrix 0O0000100a%+4+y 
LO-yr 02 ay ce ra 
Displacement Ns=|01 2 Oy? oa ry = 
interpolation 001000 —z y = 
matrix 
Geometric ma-| A= [Noi Nos Nex] 
trix 


4.4 Belarbi and Bourezane (2005-first triangular 
element) 


As mentioned previously, most of the available triangular strain-based plane 
elements have the geometry similar to the one in Table 1. In fact, many of 
these elements attempted to improve the performance of element suggested 
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Table 2: Details of the triangular element proposed by Sabir and Sfendji [72] 
Ref | Properties Geometry 
Geometry 


WN 
° 
< 


Ex (,Y) = (Ex)o + (Ex,y 


aa Strain field Ey (x,y) = (ey), + (Ene) 2 
= Yay (2, y) = Gane 
g |Strain state vec-|S= {uo Uo To (Ex), (Ey), 
A | tor Gea (Ena). Crha 
a Nodal displace- D= {Doi-1 Do; Doj-1 Do; Dop-1 Dox 
. ment vector Dg, Da }" 
4 000100y0 
“#2 | Strain interpola-;B,= |O00001002 
tion matrix 00000100 


10=y20% ay —% 


Displacement N,= 2 
Ol x Oy5- > xy 


interpolation 
matrix 

Geometric ma-| A= [Nei Nsj Nex Nei] 
trix 


by Sabir [71]. In one of these research works, Belarbi and Bourezane [9] 
proposed a new element by incorporating Poisson’s ratio in the assumed 
strain field. They suspected that unsatisfactory performance of the element 
proposed by Sabir might be due to the existence of coupling terms in the 
direct strains. Therefore, they included the Poisson’s ratio in their assumed 
strain-field (see Table 4). Like the previous strain-based triangular elements, 
the utilized strain field in this study satisfied the compatibility condition, but 
the equilibrium equations were not considered in this formulation. 
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Table 3: Details of the triangular element proposed by Tayeh [75] 
Ref | Properties Geometry 
Geometry Same as Table 1 
Z 
Ex (2,4) = (ez), a (Say) oY as (ow) t 
Gh ain Geld Ey (x,y) = (Ey), ca (€y,2) 5¢ = (Cains or 
Yay (fu) s (eeg)., + (€y,20) 5% (Syoe)Y 
—(Ezy)oa + (Ey,2) 05 
Strain state vector) S = {uo Vo To (ae eyo (tea), ae ls 
1 (Ey) (Ex,yy)ot 
&, |Nodal — displace-|Same as Table 1 
‘S| ment vector 
> z 
& 000100 y O — 
Strain interpola-|/B;= |000010 0 « -—*% 
tion matrix 000001-" ¥a+y 
10-yxr0$% ay _2 5 “4 t+ % 
Displacement Ns=|01 2 Oy g 7 e wy a ou 
interpolation 001000 -2-2 u +y% y ay 
matrix 
Geometric matrix |Same as Table 1 
4.5 Belarbi and Bourezane (2005- second triangular 
element) 


In 2005, Belarbi and Bourezane [10] performed another study and proposed 
a triangular strain-based element with the geometry similar to their previous 
work, but with a different strain field. They assumed linear variation of 
normal strain with respect to the perpendicular direction for this element, 
while the incomplete second-order field was assumed for shear strain. Once 
again, the three-node nine-degrees of freedom geometry was considered for 
this element. Further details of the element are provided in Table 5. 


4.6 Rezaiee-Pajand and Yaghoobi (2014- first triangular 


element) 


Rezaiee-Pajand and Yaghoobi [65] proposed a five-node triangular element 
with a complete linear strain field, which its geometry and details are pre- 


sented in Table 6. 


In the assumed strain field, there are nine strain states and three rigid 
body modes, which totally become twelve strain states. The selected strain 
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Table 4: Details of the first triangular element proposed by Belarbi and Bourezane [9] 
Ref | Properties Geometry 


Geometry Same as Table 1 
Ex (By) = (Ex), a (2ny) o¥ _ (€x,2),0*5% 
—(eye) 24 
& | Strain field Ey (2, y) = (Ey), + (Ey) — (Ex,0) oY 3 
a —(Ex,y)oyY 
® Yay (x,y) = (Ve) r (Sanit at (Scene) 
6 [Strain state|/S={u. vo ro (Ez), (Ey), (ey), (Ezy), 
: vector (uc). ee). ) 
I Nodal dis- | Same as Table 1 
‘3 | placement 
& | vector 
aa) 000100 y -vx = 
Strain interpo-|B,= |O00010-vy a2 = 
lation matrix 000001 0 O x+y 
. LO-y2Os oY ; v vee e 
Displacement N,= 1/01 2 Oye —S ry eH ee 
interpolation 001000 —2 y au 
matrix 
Geometric ma- | Same as Table 1 
trix 


components satisfy the compatibility requirement. In addition, Rezaiee- 
Pajand and Yaghoobi enforced equilibrium conditions in this strain field. 
They found that the necessary condition for satisfaction of the equilibrium 
criteria is that some strain states be dependent to others. Therefore, they 
selected the two dependent strain states, and as a result; ten independent 
strain states remained. In agreement with the number of independent strain 
states, the resulting element needs ten degrees of freedom. Rezaiee-Pajand 
and Yaghoobi considered a triangular element with six nodes, as depicted in 
Table 6. It is evident, four of these nodes have two translational degrees of 
freedom, while the other two only have one translational degree of freedom 
perpendicular to the corresponding side of the element. 

The displacements of mid-side nodes, which are perpendicular to the el- 
ement sides, can be connected to the displacements of the nodes in x and y 
direction, by using the following relationship: 


w=u-sin(a) +v-cos(a). (33) 


In which, a is the angle between the degree of freedom normal to the side 
and the z axis. In this case, due to the difference in the number of degrees 
of freedom at the nodes, the geometrical matrix is a bit dissimilar to the 
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Table 5: Details of the second triangular element proposed by Belarbi and Bourezane 
[10] 


Ref | Properties Geometry 
Geometry Same as Table 1 
€z (a,y) = (Ex). a (€x,y) o¥ 
Strain field yh 0) = ley) oF Cyl 


Yay e y) = (Yeu) 15 (Alay cess) 0 5 tein 
Strain state) S={uo vo To (éz), (Ey), (ev), (Ezy), 


vector Gyels (eaes)}” 
Nodal dis- | Same as Table 1 
placement 

vector 


O000100y0 O 
Strain interpo-|B,= ;|O000100a2 0O 
lation matrix 0000010027+y? 


Logos 2y =2- = 


3 
Displacement |N,= |01 2 Oy — xy = 


Belarbi and Bourezane [10] 


2 2 


interpolation 001000-2 y i 
matrix 

Geometric ma- | Same as Table 1 

trix 


previously reviewed finite elements (see Table 6). Here, the first sub-matrices 
are derived by replacing the node coordinates in the matrix N, while for 
the last two sub-matrices, a different relation was utilized (see Table 6). 


Despite its irregular degrees of freedom, the numerical evaluations showed 
that the performance of this element is better than many of the previous 
membrane elements. There are different reasons for the better performance 
of this element. First, all the components of its strain field had the complete 
term for a linear approximation. Moreover, the equilibrium equations are 
imposed on the assumed strain field, and the independent strain states are 
excluded from the strain state vector. The later property provides also the 
advantage of reducing the number of degrees of freedom in the element. These 
investigators realized that the geometry of the element and its node locations 
also had significant effects on its performance. This influence resulted in a 
new series of more accurate elements, which will be discussed in the coming 
sections [58, 59, 61, 62]. 
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Table 6: Details of the first triangular element proposed by Rezaiee-Pajand and 
Yaghoobi [65] 
Ref | Properties Geometry 
Geometry 


Rezaiee-Pajand and Yaghoobi [65] 


Strain field 


Strain state 


4 Vo To (Ex) (Ey). (Yeu) (Ex,0) 9 (Ex,y)o 


vector wtly (Eyy)ot” 
Nodal _ dis- ={Dai-1 Dai Dej-1 Doj Deor-1 Dor Doai-1 Dai-1 
placement i sD? 
vector 
Strain in- | B, 
terpolation 000100 x y 0 0 
matrix = |000010 0 0 x y 
Q2G+A)y _ d dN. @G4A)« 
Govgor1 Oey ty Ay Cer 
Displacement N, 
interpola- _ |10-ya0 ¥ ha ACNE xy —¥ (3g + 2) , 2 ; 
tion matrix Ol a Oys 0 —2? Gg + 3) ry x _ OGtr)e 
Geometric | A= [N5i Noj Nex Noi Nem Non] 
matrix Nsg 
= a» (a) cos(a) agcos(a) — ygsin (a) 
agsin(a) —ygcos (a) ee _— yasin(a) | 
2 2 
(3 = ame ) sin (a) 
(xeyg)sin(a) — (xg? (se + 3) cos (a) 
(xgyg) cos (a) — (ys? ae + 3)) sin (a) 
2 2 
uy — Cor 22") cos (a) =m,n 
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4.7 Rezaiee-Pajand and Yaghoobi (2014- second 
triangular element) 


In another study, Rezaiee-Pajand and Yaghoobi [66] utilized the complete 
linear strain field of the previous study, but with a different element geom- 
etry. In this work, they proposed a seven-node triangular element, which is 
depicted in Table 7. Six of the ten required degrees of freedom were allo- 
cated at the vertex nodes, which had two translational degrees of freedom 
each. Three mid-side nodes had only one translational degrees of freedom 
perpendicular to their sides, and the last degree of freedom was a drilling one 
at the center node. The first three sub-matrices for nodes i, 7, and k of the 
geometric matrix, A, are derived by replacing the coordinates of these nodes 
in the displacement interpolation matrix. For the three mid-side nodes, the 
respective sub-matrix is computed by using the equation given in Table 6 
for N,g. Finally, the last sub-matrix, Ns,, is derived from the next equality 
presented in the last row of Table 7, after the geometric matrix. The authors 
developed this formulation for geometrical nonlinear analysis of plane struc- 
tures. For this purpose, they took advantage of the co-rotational formulation 
(54, 3]. For this purpose, a local coordinate system, which its origin was 
located at the center of the element, was considered. This coordinate system 
translates and rotates with the element. Then, the element was formulated 
in this new system. Since the nonlinear analysis was not in the scope of this 
study, for further information about the co-rotational formulation, one can 
refer to references [54, 66]. 


4.8 Rebiai (2018) 


In a more recent attempt to propose three-node nine-degree of freedom trian- 
gular element, Rebiai [48] suggested a new strain-based element with incom- 
plete second-order strain field. This strain field satisfies the compatibility 
condition, but the equilibrium conditions were not fulfilled within the ele- 
ment. In fact, Rebiai did not intend to impose the equilibrium equations, 
(Please see Table 8). 


4.9 Rezaiee-Pajand, Gharaei-Moghaddam, and 
Ramezani (2019- first triangular element) 


In one of the most-recent studies, Rezaiee-Pajand, Gharaei-Moghaddam, and 
Ramezani [58] suggested new triangular elements. They utilized the complete 
linear strain field, which was the same as the previous studies by Rezaiee- 
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Table 7: Details of the second triangular element proposed by Rezaiee-Pajand and 


Yaghoobi [66] 


Ref | Properties 


Geometry 


Geometry 


—______»> ' 


A= 
Xp [ 01000 


Ss Strain field Same as Table 6 

‘= | Strain state | Same as Table 6 

9 vector 

an Nodal dis- D= {Doj-1 Do; Doj-1 Do; Dorp Doz D, 

Ss placement Dm Dn Det 

z vector 

® Strain in- | Same as Table 6 

zg terpolation 

= matrix 

7 Displacement| Same as Table 6 

5 interpola- 

N tion matrix 

| Geometric [Nei Noy Nox Noi Nem Non Nop] 
matrix 


(2G+A) yp 
“3G 
(2G+A)tp  (2G+A) yp (2G+A)ap 
2G 2G IG 


Pajand and Yaghoobi [65, 66]. Moreover, they also imposed the equilibrium 
condition to specify the dependent strain states. The differences between 
these new elements and previous ones are in the geometry of the elements. 
In their first element, the authors assumed a five-node triangular element, 
with three vertex and two mid-side nodes (see Table 9). Each node has two 
translational degrees of freedom. 


As mentioned previously, it has been proved through extensive investi- 
gations, that configuration of the nodes and distribution of the degrees of 
freedom have been influential in the performance of the resulting elements 
(65, 66]. Numerical studies showed that the last presented element provided 
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Table 8: Details of the triangular element proposed by Rebiai [48] 
Ref | Properties Geometry 
Geometry Same as Table 1 
éx (2,y) = (Ex) + (Ex,y) oy + (Ex,yy) oY” 
Strain field éy (v,y) = (Ey), + (Ey,a) 2 + (Ex,yy) 2? 
Yey (@,Y) = Yey)o + 2(Ea,yy)o%” + 2(€a,yy) oY" + 4(Ex,yy)oLY 


‘oO 
=H oO 
vq [Strain state]S={uo vo ro (Ex) (Evlo Mevlo (Eaulo (eux)o 
‘3 | vector (Ex,yy)}7 
wo Nodal _ dis- | Same as Table 1 
placement 
vector 


00010070 ¥ 
Strain inter-|B,;=|00001002_ 2? 


polation ma- [00000100 2(«+y)? 
trix 
y y> 2? 2 

10-yrxO 5 zy 5 3, 7 78 
Displacement Ns=]01 2 0y Z_f gy a + yx? 
interpolation 001000 -x y x-y 
matrix Zz 
Geometric Same as Table 1 
matrix 


very accurate responses, especially for the elements under the bending loads. 
Something that might be questionable about the geometry of this element 
was the arbitrariness in the selection of the sides with mid nodes. To show 
that this selection did not have considerable effect on the performance of the 
element, the authors solved a numerical example with different locations of 
the mid-side nodes, and the results were identical. Moreover, the main idea 
behind this selection was to produce a powerful finite element for transitional 
purposes. To demonstrate this fact, Rezaiee-Pajand, Gharaei-Moghaddam, 
and Ramezani [58] investigated the performance of this element as a transi- 
tional element in numerical example and compared its behavior with standard 
displacement-based transitional elements. The results showed superiority of 
the mentioned element. 


4.10 Rezaiee-Pajand, Gharaei-Moghaddam, and 
Ramezani (2019- second triangular element) 


After presenting the previous element, Rezaiee-Pajand, Gharaei-Moghaddam, 
and Ramezani [58] considered another configuration for the complete linear 
strain field. In this element, the authors considered the well-known three node 
nine-degree of freedom triangular element and added an internal node with 
one translational degree of freedom in an arbitrary direction. This element 
geometry is demonstrated in Table 10. By using (33), the displacement of the 
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Table 9: Details of the first triangular element proposed by Rezaiee-Pajand, Gharaei- 
Moghaddam, and Ramezani [58] 
Ref | Properties Geometry 
Geometry 


Strain field | Same as Table 6 
Strain state | Same as Table 6 
vector 
Nodal dis- | D= {Doi-1 Dy; Doj-1 Do; Dox-1 Dox Dai-1 Dai 
placement Dom-1 Dom}? 


Rezaiee-Pajand and Gharaei-Moghaddam [58] 


vector 
000100 x y 0 0 
Strain in- | B, = |000010 0 0 x y 
terpolation 000001 Qa aly dar » Gee 
matrix " 7 
Displacement N, 
interpola- _ [10-yard $ a QE yy xy -y? (sj +4) 0 
tion matrix Ol a Oys 0 —a? (A, +4) xy ve _ ae 
Geometric | A= [Nei Ns; Nex Noi Nem] 
matrix 


internal node can be connected to the displacements in x and y directions. 
The internal node of the element can be removed by the static condensation 
approach, and therefore, the element becomes a three-node nine-degrees of 
freedom triangular element, which is a common element in general finite el- 
ement programs. This element provided surprisingly accurate results under 
different types of loading and especially shear loading, in which many of the 
available elements failed to provide the exact response in the case of distorted 
mesh [58]. In addition to the fast convergence and high accuracy, this ele- 
ment was highly insensitive to the mesh distortion. Numerical examinations 
demonstrated that the direction of internal degree of freedom did not have 
any noticeable effect on the accuracy and performance of the element [58]. 
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Table 10: Details of the second triangular element proposed by Rezaiee-Pajand, 
Gharaei-Moghaddam, and Ramezani [58] 

Ref | Properties Geometry 

Geometry 


2 Strain field | Same as Table 6 
- Strain state | Same as Table 6 
ES vector 
Z | Nodal dis- | D={Dy;-1 Dai Daj-1 Daj Dop-1 Dor Dai-1 Dai 
= | placement. Dom—1 Dom }" 
s vector 
3 000100 x y 0 0 
5 Strain in- | B, = |000010 0 0 x y 
oy terpolation 000001 OGtNy Aw Ay GG tA)s 
Z matrix 
8 Displacement Ng 
Bie acel_ ee e  ) : e 
| eee ee el dee ee eee 
g 00 a 1 00 0 3G = 2G “2 2G ue 3G “2 
‘3 B=i,j,k . 
& | Geometric | A= [Nei Nsj Nex Noi] 
matrix 
Ns, = |sin(a) cos(a) xcos(a) —ysin(a) «sin(a)  ycos (a) 
* 2 
sitesia) ysin{a) ar GGr Nu ) sin (a) 
(aim) sin (a) — (21? (34 + §)) co (a) 
(xim) cos (a) — (wi? (34 + §)) sin (a) 
(8a) 


4.11 Rezaiee-Pajand, Gharaei-Moghaddam, and 
Ramezani (2019- third triangular element) 


In another study, Rezaiee-Pajand, Gharaei-Moghaddam, and Ramezani [59] 
formulated a higher-order strain-based assuming complete second-order nor- 
mal strains and linear shear strains (see Table 11). After imposing the equi- 
librium and compatibility conditions on the assumed strain field, seven de- 
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pendent strain states were excluded and eleven independent strain states 


remained. 
Table 11: Details of the third triangular element proposed by Rezaiee-Pajand, 
Gharaei-Moghaddam, and Ramezani [59] 

Ref | Properties | Geometry 

Geometry 
f oy. 
2 
Ex (x,y) (Ex)o t (Ex,0) 5 t (x,y) oY (Ea,00)g (=) 
2 
+(Ez,2y), (xy) + (Ex,u)o (4) 
= “ 7 a 
ZF | Strain field | 9 cy (x,y) = (ey), + (Cywlot + Cv )ov + (Cxaz), (5) 
* at 

S +(€y,2y), (wy) + (Evy) o (4) 

i Yay (x, y) = (Yey)o + (Yey,2) ot + (Yeyy)oY 

& | Strainstate | S={uo vo ro (Ex), (Ev)o Ory)o (Exe)o (Exy)o (Eve) 

= vector (Eyy)o (Exyy) ot” 

b Nodal  dis- D= {Doi-1 Do; Doj-1 Do; Dox Dox D, Dm Dy 

oO 

a placement Dop-1  D2p} 

2 vector 

Strain in- | B, 

. 2 Od 
ee 000100 « y 0 0 w+ AGCEN) 
matrix = x? BY. 

=!000010 0 0 « asey = = ak 
2G+A)a d 2G+A)x 
000001 —GGtAw _Ag Ay —OGt 0 
Displacement N, 
interpola- — f10-yro ee - CGAY gy —y?( GA) 0 (acicx) + ay 
tion matrix Ol x Oys 0 2 G+) : zy ye (2G ae a aocx) yx? 
Geometric A= [Noi Ns; Nsk Tst Tsm Tsn Neo| 
matrix T.- 


(2G+4+A)y — (2G+A)# (2G+A) 
[001000 cA ae tal 


Q2G+4A)x 
IG 


The selected geometry of the element is demonstrated in Table 11. As it 
can be seen, the element had seven nodes and eleven degrees of freedom in 
agreement with the independent strain states. 
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4.12 Rezaiee-Pajand, Gharaei-Moghaddam, and 
Ramezani (2020- fourth triangular element) 


Rezaiee-Pajand, Gharaei-Moghaddam, and Ramezani [61] utilized the as- 
sumed strain of their previous study [59] to formulate another higher-order 
strain-based element. The difference of this element with the previous one is 
the geometry of the element, which is demonstrated in Table 12. 


Table 12: Details of the fourth triangular element proposed by Rezaiee-Pajand, 
Gharaci-Moghaddam, and Ramezani [61] 
Ref | Properties 


Geometry 


Strain Same as Table 11 
field 


Strain Same as Table 11 
state 


vector 
Nodal dis- | D= {D3i_2 D3-1 D3: D3j_2 D331 D3; D3x—2 
placement | D3,-1 D3, Do; D2}7 

vector 
Strain Same as Table 11 
interpo- 
lation 
matrix 
Displacement Same as Table 11 
interpo- 
lation 
matrix 
Geometric | Ag = [Nei Tsi Noj Taj Nox Ts: Nor Nom Non Neo] 
matrix 


Rezaiee-Pajand et al. [61] 
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4.13 Rezaiee-Pajand, Ramezani, and Gharaei- 
Moghaddam (2020- fifth triangular element) 


Since it was expected that using higher-order strain fields improves the accu- 
racy of the elements, Rezaiee-Pajand, Ramezani, and Gharaei-Moghaddam 
[62] formulated a seven-node triangular element using complete second-order 
strain field (see Table 13). In this element formulation, the compatibility and 
equilibrium conditions were again satisfied. 


Table 13: Details of the fifth triangular element proposed by Rezaiee-Pajand, 
Ramezani, and Gharaei-Moghaddam [62] 

Ref | Properties 
Geometry 


Ex (t,y) = (Ex), + (€x,n) 9% + (Ex,y)oy + (Exc), (2) 


2 
s y 
+(Ex,cy), (xy) ae Eawy)s (4) 
— a 
= ey (2,9) = (Ey), + Eyc)ot + (Cv )o¥ + (Ever), (F 
ag Strain field y 
~ +(Ey,cy), (wy) + (Ey,yy)o (4) 
=| 
a _ r 
Yay (x, y) = (Yey)o + (Yey,x ot + (Yeyy)oY + (Yey,02)y (=) 
2 
ne ¥ 
5 +(Yey,2y)o (ay) + (Yey.vu)o (4) 
& Nodal dis- | D={Dg-1 Dx; Daj-1 Do; Dor-1 Dor Dai-1 Dai 
P T 
b placement Dom—1 Dom Dan—1 Don Doo-1 Doo} 
3 vector 
s Strain in- | B, 
P Pa eee Sioa 
al terpolation 000100 & y 0 0 HCH’ _ x ay alan 0 
‘ = (Gir) 22 Gy ; 
matrix te ee ee y eae age FEY Ss 
000001 ~% ~g2 -ay -¢ vy =e 3G tye - 3G 


interpola- lo-yroy se —-2Grwh gy (GA) 

tion matrix 
- T 

Geometric | Ag = [Noi Nj Nsz Nsi Nem Nen Noo] 


matrix 


Displacement NV, 


yo (2G+4)y" Gtr Gtr 
up — ECe -(S) WG) 


ny? _ (2G4+A)2" 
2 re 


Ole Oye 0 —2?(GtA) ay = CGH 
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5 Quadrilateral membrane elements 


The existing quadrilateral strain-based membrane elements are reviewed in 
this section. In contrast with the triangular elements, which are mostly 
developed by the assumption of linear or second-order strain-field, higher- 
order strain fields were utilized in the formulation of quadrilateral elements. 
In the quadrilateral elements, there are two common configurations. The 
first one is a four-node twelve-degrees of freedom quadrilateral element, and 
the second common geometry is a five-node ten-degrees of freedom element. 
More details are provided in the following. 


5.1 Sabir and Sfendji (1995) 


Like the triangular elements, the first quadrilateral element is proposed by 
Sabir [72]. In 1995, Sabir and Sfendji suggested a rectangular strain-based 
finite element with linear strain field. 


Once again, they assumed strain components satisfy compatibility con- 
dition, but the equilibrium equations are not fulfilled. Due to the existence 
of ten strain states, ten degrees of freedom are required. Therefore, Sabir 
and Sfendji considered the five-node rectangular element depicted in Table 
15. Each node of this element has two translational degrees of freedom. 


5.2 Tayeh (2003) 


Another rectangular element was proposed by Tayeh [75]. He utilized an 
incomplete fourth-order approximation for normal strains and incomplete 
second-order assumption for the shear strain (see Table 15). 


In agreement with these twelve strain states, a four-node rectangular ele- 
ment with three degrees of freedom at each node was considered. It should be 
noted that this element and the previous one proposed by Sabir and Sfendji 
had rectangular geometry, according to the original articles [72, 75], but 
based on the basics of strain- formulation, the effect of the geometry only 
entered in geometric matrix, A. So, by using the nodes of a general quadri- 
lateral shape in construction of this matrix; the mentioned elements can have 
a general quadrilateral geometry. 
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Table 14: Details of the Rectangular element proposed by Sabir and Sfendji [72] 
Ref 


Properties 
Geometry 


tion matrix 


Los, 
a Strain field | Same as Table 2 
v= | Strainstate | S={uo vo To (Ex), (Ey)o (Yay)o (Exylo (Ey,2)o 
zg vector (Yey,2)o Gara 
a Nodal  dis- D= {Doi-1 Do; Doj-1 Do; Dop_4 Do, Do 
vv placement Dy Dam—1 Dom }" 
S vector 
2 000100y000 
a Strain in- | B, = |0000100200 

terpolation O00001L100ay 

matrix 

Sf 2 
Displacement N, = 1O=yens ae > is 2 7 
interpola- Ole Oys—-5 ary “0 


Geometric 
matrix 


A= [Nei Ns; Nsz Noi Nem] 


5.3 Belarbi and Maalem (2005) 


In 2005, Belarbi and Maalem [12] suggested an improved strain-based rect- 
angular element by the assumption of linear normal strains and constant 
shear strain. The element geometry and distribution of nodes and degrees of 
freedom are the one proposed by Sabir and Sfendji. 
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Table 15: Details of the Rectangular element proposed by Tayeh [75] 


Ref | Properties 
Geometry 
AOn-s 
Ds, 
_ 
A 
>* 
Ex (t,Y) = (Ex), + (Ex,y) oY 4 (Ex,yy) oY? + (Ex,cyyy)92y° 
Strain eld Ey (2, y) = (Ey), + (Ey,0) 9% — (Ex,yy) 2? — (Exam) 
= Jey (%,Y) = Maule + O2v.2)ot + Yeu oy — (Eru)or 
1 
a (Ey, °F 
g [Strainstate |S={uo % To. (xlo Evo Cav), Evade Evade 
& vector (AEey«)y Oevy)o (Exyy)o (Ex.2yyy)o 
Nodal  dis- D= {D3i-2 D3i-1 D3; D3;-2 D3j-1 D3; D3p_2 
placement | D3x-1 D3z D3i-2 Dsi-1 Dsi}" 
vector 
Strain in- | B, 
terpolation 000100 y 0 00 y? ay? 
matrix —=|000010 0 «@ 00a? —xy 
2 
000001-S-¥ zy 0 0 
Displacement Nx 
A 3 42 2 2243 
ete 10-yar0 4 sy ae 0 yay? 2 : 
10n Matrix = Olgz O¥s-5- 5 wy £ 0 xy en 
001000 -S-2 Y+y § -¥ -2ay 354 
Geometric A= [Nei Ns; Nex Nail 
matrix 


5.4 Rezaiee-Pajand and Yaghoobi (2012- first 
quadrilateral element) 


The first strain-based generalized quadrilateral plane element, which its strain 
field satisfies both compatibility and equilibrium conditions, is proposed by 
Rezaiee-Pajand and Yaghoobi [63]. They [65, 66] took advantage of the lin- 
ear strain field that they had used for development of triangular elements 
(see Table 6 ). After imposing the equilibrium criteria and determination 
of dependent strain states, ten independent strain states remain for the re- 
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Table 16: Details of the Rectangular element proposed by Belarbi and Maalem [12] 
Ref | Properties 
Geometry Same as Table 14 


Ex (2, y) = (Ex) . + (Ex,y) oy + (Ex,x) 5% 


ay Strain field Ey (2, y) = (Ey), + (Ey,2) 8 + (Evy) oY 
g Yay (z,y) = (Yey)o 
< Strain state]S={uo vo To (Ex), (Ey), (eyo (Exy)y (Ex,x), 
s vector (Eye) (Ey.y)ot™ 
a Nodal — dis- | Same as Table 14 
a placement 
3 vector 
a O00100yzxz00 
a Strain inter- |B; = f 00010002 
polation ma- 0000010000 
trix 
y oy 
Displacement | NV; = ; an oe ie 2 2 | 
interpolation Ola 0y5-5 0 ay % 
matrix 
Geometric Same as Table 14 
matrix 


quired ten degrees of freedom. Rezaiee-Pajand and Yaghoobi considered the 
generalized five-node quadrilateral element, which is depicted in Table 17. 

In another study, Rezaiee-Pajand and Yaghoobi [64] investigated the per- 
formance of two special rectangular variants of the generalized quadrilateral 
element. In the first element, the fifth node was located in its center, while 
this node was moved to the middle of one side in the next element. In this 
study, they showed that these strain-based elements were less sensitive to 
mesh distortion in comparison with their displacement-based counterparts. 
Moreover, they displayed that strain-based elements were completely free 
from shear parasitic errors. 


5.5 Rebiai and Belounar (2013- first quadrilateral 
element) 


A four-node rectangular strain-based element with incomplete fourth-order 
normal strains was proposed by Rebiai and Belounar [49]. The authors as- 
sumed linear shear strains for this element. This strain field, had twelve 
independent strain states. Regarding this field, it seems that there was no 
clear and rational basis behind this selection. Therefore, this field was prob- 
ably the result of the trial-and-error process to attain the best performance 
of the resulting element. Rebiai and Belounar developed this element for ax- 
isymmetric, as well as, elastoplastic analysis of structures. For the nonlinear 
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Table 17: Details of the first quadrilateral element proposed by Rezaiee-Pajand and 
Yaghoobi [63] 

Ref | Properties Geometry 
Geometry 


Strain field | Same as Table 6 
Strain state | Same as Table 6 
vector 
Nodal dis- | D= {Doi-1 Do; Doj-1 Do; Dop-1 Dox =Dai-1 
placement Dy Doam—1 Dom" 

vector 
Strain in- | Same as Table 6 
terpolation 
matrix 
Displacement Same as Table 6 
interpola- 
tion matrix 
Geometric A={Ns Nei Nex Nat Nowak 
matrix 


Rezaiee-Pajand and Yaghoobi [63] 


analysis, they took advantage of the Mohr—Coulomb yield criteria and the 
initial stress and strain methods [49]. They showed that the strain-based 
elements provided reliable and accurate results in nonlinear problems. 


5.6 Rebiai and Belounar (2014- second quadrilateral 
element) 


In another research work, Rebiai and Belounar suggested another variant of 
their previous element [50]. In this new element, they considered the strain 
field of their previous study [50], but added new linear term to the shear 
strain and changed the dependent term of the normal strains. 
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Table 18: Details of the first quadrilateral element proposed by Rebiai and Belounar 


[49] 
Ref Properties 
Geometry i 
Ex (t,y) = (Ex), ote (Ex,y) oY + (Ex,x) 9% + (Ex,uy)o¥* 
+2(Ex,cyyy)otY* 
Ey (2, y) = (Ey), + (€x,2) 52 
gz Strain field +2 (Ey,y) oY = (En,yy) 9%" 
a —2(Ex, yyy )oYE? 
: nay (9) = 2 ey)o + 2Enw)o@ + 2Ene)o¥ +2 (Ew )oy 
Ee +2(Ex,yy) oY = 2(Ex,yy) © +2 (eye) 
=~ | Strainstate | S={u, vo To (Ex), (Ey), (Mey), (Exy), (Exc), 
a vector (Ey, ys (Yey,x)o (Exau)o (Ex,evyy)ot 
‘ge | Nodal dis- |] D={D3i-2 Ds3i-1 D3: Dsj-2 Dsj-1 Ds3j Dsr—2 
2 placement | Dsp-1 D3 Dsi-2 Dsi-1  D3i$* 
vector 
000100y 200 ¥ Qry? 
Strain ine | B,=]}0000100 2 2y 0 —-a? —2r3y 
terpolation 00000 2 2x 2y 2y 2a 2Qy-2r 0 
matrix 
Displacement N, 
interpola- 10-yx0yay wey y2 0 aye+y? a2y3 
tion matrix |= |9 1 7 QO yx . ry yr a —x2y—2? ay? 
0010000 0 y © —a@—y—2xy —3x7y? 
Geometric A= [Nei Ns; Nex Neil 
matrix 


Similar to their previous study, Rebiai and Belounar utilized this element 
for both linear and materially nonlinear analysis. In their numerical evalua- 
tions, they included different yield criteria, such as, Tresca, Von Mises and 
Mohr-Coulomb. Once more, the attained results demonstrated the superi- 
ority of strain formulation in comparison with classical displacement-based 


elements. 
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Table 19: Details of the second quadrilateral element proposed by Rebiai and Belounar 

[50] 

Ref Properties 

Geometry Same as Table 18 

Ex (2, Y) = (€x)o + (Ex,y) oY + (Ey) ¢ + (Ex,yy) oY? 
+2(Ex,2yyy) ty? 


Ey (2, y) = (Ey), +r (Ey,0) 5© +r (Ey,y oY = (aig) 


S | Strain field é 

— (ery AY ad 

A Yay (2, Y) = 2 (Yay) + [2(Yay)o + 2(Ey,c)5 + (Evy) ol¥ 

8 +2[(Coino + Evo + Caudal? 

M | Strainstate | S={uo Vo To (Ex), (Ey), (Wey), (Exy), (Ey,0)5 
= | vector ei), Gada): Gea Eman 

q ¥Ulo LY,t)o z,yylo Z,XYyylo 

s Nodal dis- | Same as Table 18 

s placement 

oe vector 


0001s 0 yoOoOO y day? 
Strain ine | B,=]00001 O O wy O —a? —-22%y 


terpolation 0000 2a 2y+2 2x 2yy2x 0 0 

matrix 

Displacement N, 

] le 12 2 2,2 

interpola- 10-ye = ¥? ty wy a z 0 ay? ay? ] 

tion matrix | — Ol xr 0x2 +y 2 a zy v x2 xy wy 
0010 «2 -y 0 0 -§ & —2ay —3a*y 

Geometric Same as Table 18 

matrix 


5.7 Rebiai, Saidani, and Bahloul (2015- Third 
quadrilateral element) 


In the following of his previous researches, Rebiai, Saidani, and Bahloul [51] 
suggested another strain-based quadrilateral element for linear dynamic anal- 
ysis of the plane problems. They utilized the general form of the strain field, 
which was used in their previous studies, but with slight modifications (see 
Table 20). The element geometry and its strain state and nodal displace- 
ment vectors were similar to the previous element proposed by Rebiai and 
Belounar [50]. The difference between these two membrane elements was in 
their strain and displacement interpolation matrices. Application of this new 
element for dynamic analysis of the plane problems proved acceptable accu- 
racy of the assumed strain method for the dynamic problems. They utilized 
the lumped mass matrix for the element. 
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Table 20: Details of the quadrilateral element proposed by Rebiai, Saidani, and Bahloul 
[51] 
Ref Properties 
Geometry | Same as Table 18 
Ex (%,y) = (Ex), + (Ey), ot (Ex,y) oY oT (Ey) © a (Cay) sur 
+2(Ex,cyyy)92y° 
— | Strain Ey (L,Y) = (Ey), + (Ey) 9% + (Ey,y) oY — (Ex,yy)o2" 
= 3 
2, | field —2(Ex,2yyy) YX 
Z Jey (£,Y) = 2 (Yay), + [2Oyan ts + 2(€y,2), + (ey), y 
aS +2 [aa ey Orie). | x 
‘Ss | Strain S={uo Uo To (Ex), (Ey), (ey), (Exy), (Eve), 
a T 
oe state (Eyyo Oeyelo (Exyy)o (Ex,cyyy)ot 
vector 
Nodal dis- | Same as Table 18 
placement 
vector 
000laz+l 0 y O00 xy 2x 
Strain B,=|0000 1 0 0 xy 0 —a? —223y 
interpo- O000 2 2Z+222xayy2xr 0 0 
lation 
matrix 
10-yx = +a y?4 y ry v 5 0 ay? ay 
Displacement Ns = |01 4 O0a?+y «& tay a? —a?y —a3y? 
interpo- 0010 « -y 0 0 -% & —22y —327y? 
lation 
matrix 
Geometric | Same as Table 18 
matrix 


5.8 Rezaiee-Pajand and Yaghoobi (2015- fourth 
quadrilateral element) 


In most of these formulations, there is no rational basis for the selection of 
the assumed strain field. However, Rezaiee-Pajand and his colleagues pro- 
posed application of the concept of Taylor expansion for this purpose. In an 
attempt to develop second order strain-based elements, Rezaiee-Pajand and 
Yaghoobi [67] suggested the complete second-order strain field (see Table 21). 
Excluding these dependent coefficients from the total strain states, fourteen 
independent strain states were remained. Accordingly, in this element, the 
four vertex nodes had three degrees of freedom, namely two displacement and 
one in-plane rotation, and the internal node had two translational degrees of 
freedom. 
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5.9 Rezaiee-Pajand and Yaghoobi (2015- fifth 
quadrilateral element) 


Rezaiee-Pajand and Yaghoobi also proposed another second-order element, in 
which the assumed strain field was the same as the one presented in Table 21. 
In addition, the equilibrium criteria were imposed only on the linear terms 
of the strain components [67]. Therefore, only two strain states were charac- 
terized as dependent ones, and eighteen independent strain states remained. 
To generate an element with eighteen degrees of freedom, they considered a 
nine-node generalized quadrilateral element, which is demonstrated in Table 
22. Each node of this element had two translational degrees of freedom. 


5.10 Hamadi, Ayoub, and Maalem (2016) 


In 2016, Hamadi, Ayoub, and Maalem [31] independently proposed a new 
quadrilateral finite element. This element is indeed the same as the element 
previously proposed by Rezaiee-Pajand and Yaghoobi [63], restricted to rect- 
angular shapes. 


5.11 Rezaiee-Pajand and Yaghoobi (2018- sixth 
quadrilateral elements) 


In order to analyze geometrically nonlinear plane structures, Rezaiee-Pajand 
and Yaghoobi [63] modified their five-node quadrilateral element by the co- 
rotational approach [70]. Their findings showed that the strain-based for- 
mulation can provide accurate results for geometrically nonlinear analysis 
of structures. Besides, it did not show any sensitivity to the aspect ratio 
and mesh distortion. To sum up, the summary of the reviewed elements is 
presented in Table 23. 


6 Other types of strain-based elements 


Since the focus of the present study is on the membrane elements, only the 
available strain-based elements were reviewed in the previous sections. How- 
ever, the advantages of strain formulation persuade researchers to have taken 
advantage of this approach in the development of finite elements for other 
types of structures. In order to provide a brief introduction to the appli- 
cation of assumed strain technique for the proposition of the different finite 
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Table 21: Details of the fourth quadrilateral element proposed by Rezaiee-Pajand 
and Yaghoobi [67] 


Ref Properties 
Geometry 


Ex (x,y) _ (Ex), Tr Ex,x) 5 t+ (Ex, oY a (€x,22) 55 
y 


( 
+(Ex,ry) ey + (Ex,yy) oF ; 

t (Ey, a)ot Tv (Cxuo¥ ale (€y,22) o> 
+(Ey,2y) ty + (Ey,yy)or 5 
Yay (x,y) = Vay) a (Yey,) 5% a (Yeu,y)o yt Geyer) oF 

+(Ex,2y + Ey,20) CY + Geni 
Strain state] S={uo Uo To (Ex), (Ey), (May), (Exc), (Ezy), 
vector (Ey,2)o (Ey.v)o (Ex,2y)o (Ex,yy)o (Ey,cy)o (Eye) o}7 
Nodal  dis- D= {D3;_2 D314 D3j D3;~2 D35-1 D3; D3p_-9 
placement | D3.-1 D3~ Dsi-2 Dsi-1 D3: Dam—1 = Dam}" 
vector 
Strain in- | B, 
tepolation | fess gt tt Sa 
matrix 000001 —26py Ae Ay 2g — 2g? — 2a? 
Displacement JN, 
interpola- [ie ie a 
tion matrix 
Geometric | A= [Noi Toi Nsj To; Nox Tox Nor Tot x] 


matrix T. = [001000 2G4Ay —2 — Wie y Ph 2A y — 25 Aw — 2G (y? — 0?) — 9 2A (a? — y?) 


Strain field 


ei" e 


Rezaiee-Pajand and Yaghoobi [67] 


G. yn = See 
Bev Damme Gey | —aests 
Cie" Gah ve 
€ ne xy? z 


elements, a short review of the other works is presented in this section. Need- 
less to say, for more details, it is necessary to refer to the original references, 
which are cited in the following paragraphs. 


In the case of plate bending analysis, Belounar and Guenfoud [15] pro- 
posed a four-node rectangular plate element by assumption of linear curva- 
ture and second-order shear strains. The assumed strain field for the element 
only satisfies the compatibility condition. They showed that this element is 
more efficient than the corresponding displacement-based element. In 2014, 
Hamadi et al. [30] developed a rectangular element for plate bending anal- 
ysis, based on the Kirchhoff theory, and compared its performance with the 
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Table 22: Details of the fifth element proposed by Rezaiee-Pajand and Yaghoobi [67 
Ref Properties 


Geometry 


Strain field Same as Table 21 


= 
©. [Strain state[S={uo vo ro (Ex), (ey), (ayo (Ex,2)5 
ras vector (Ex,y)o (Ey,)o (Ey,y)o (Ex,22)o (Ex,2y)o T 
8 (Exyy)o (Ey,wx)y (Ey,cy)o (Eyyyo Oey,2x) (Oay,yy)ot 
g Nodal dis- | D= {Dai-1 Do; Doj-1 Do; Dop—1 
a placement Dox Dai-1 Dat Dam-1 Dem Den-1 Dean 
q vector Dop-1 Dap Doq-1 Doq Dop—1 Dop}F 
- | Strain inter- |B, 
& | polation ma- eee x y 0 0 =ay* 0000 a 
= | trix =!000010 0 0 2 000a¥ 0 0 
2 42 
$ looo001 2G Ay — Aa —Ay —2S2 0 0 ay ey 0 0 & “| 
= Displacement | NV; 
2 interpolation _|10-y20 y = 2G A y2 ay 0— GHA y2 
matrix = x G+aA 2 ‘ 
Olav Oys 0 3G xy 
3 2 3 3 
0 oa a <8 5 0 0% | 
a 3 2 2 3 3 
y 2G+A 2 x ay yu y? o 
2 ag tO g 9 3 D oe 9 4 
Geometric |A=[Ns; Ns; Ns, Ns, Nem Non Nep Neg Ns-| 
matrix 


displacement-based formulation. They showed the superiority of strain-based 
formulation in removing the shear locking problem. Another thin plate el- 
ement, based on assumed strain approach, is proposed by Abderrahmani, 
Maalam, and Hamadi [1]. In this study, they utilized higher-order strain 
field in comparison with the previous elements. In another study, Abder- 
rahmani et al. [2] formulated a new strain-based sector element for linear 
analysis of circular thin plates in 2017. Many years before this study, an- 
other sector strain-based element was proposed by Belarbi and Charif [11]. 
In 2018, the first triangular plate element, based on assumed strain approach, 
is proposed by Belounar, Benmebarek, and Belounar [14]. In this study, the 
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Numbfr Numbfr 
No. | Element Geometry of of Assumed strain field Drilling Other 
: features 
nodey _ dof 

Normal Shear Optimal 

strain strain criteria 
1 [Sabir [71] Triangular [3 9 Complete | Complete Jompatibility| yes [- 

linear linear 

2 [Sabir and Sfendji [Triangular [4 8 Tncomplete | constant Compatibility| no [Can be used as 
(72) linear transitional ele- 

ment 

3 | Tayeh [75] Triangular [3 9 Incomplete [Incomplete | Compatibility! yes [- 

second-order | second-order 

T |Bdabi and | Triangular [3 9 Complete [Incomplete | Compatibility| yes | Poisson’s ratio is 
Bourezane [9] linear linear included in the 

element 

5 | Belarbi and | Triangular | 3 9 Tucomplete [Incomplete | Compatibility) yes | - 

Bourezane [10] linear second-order 

6 | Revaiee-Pajand [Triangular [6 10 Complete | Complete | Compatibility[no [Can be used as 
and —-Yaghoobi linear linear Equilibrium transitional ele- 
[65] ment 

7 | Rezaice-Pajand [Triangular [7 10 Complete [Complete | Compatibility| yes |= 
and Yaghoobi linear linear Equilibrium 
(66] 

8 | Rebiai (45) Triangular [3 9 Incomplete [Complete | Compatibility[ yes __[- 

second-order | second-order 

9 | Resaiee-Pajand, [Triangular [5 10 Complete [Complete | Compatibility[no [Can be used as 
Gharaci- linear linear Equilibrium transitional ele- 
Moghaddam, ment 
and  Ramezani 

10 Triangular [7 10 Complete [Complete | Compatibility| yes |= 

linear linear Equilibrium 
Moghaddam, 
and  Ramezani 
58) 

TT | Rezaice-Pajand, [Triangular |7 is Complete [Complete | Compatibility : 
Gharaci- second-order | linear Equilibrium 
Moghaddam, 
and  Ramezani 
59) 

12 [Rezaice-Pajand, [Triangular [4 iat Complete [Complete | Compatibility] yes |= 
Gharaei- second-order | linear Equilibrium 
Moghaddam, 
and  Ramezani 
(61. 

13 | Rezaice-Pajand, [Triangular |7 W Complete [Complete | Compatibility| yes |= 
Ramezani, second-order | second-order | Equilibrium 
and Gharaei- 

Moghaddam 
[62] 

14 | Sabir and Sfendji | Rectangular [5 10 Complete Complete = no : 
72 linear linear 

15 | Tayeh [75] Rectangular [4 12 [Incomplete [Incomplete [Compatibility] yes |= 

fourth-order _| second-order 

16 |Belabi and | Rectangular [5 10 Complete [Constant | Compatibility[no |= 
Maalem [12] linear 

17 | Rezaice-Pajand | Quadrilateral |5 10 Complete [Complete | Compatibility[no |= 
and Yaghoobi linear linear Equilibrium 
(63] 

18 | Rezaice-Pajand [Rectangular [5 10 Complete [Complete | Compatibility[no |= 
and Yaghoobi linear linear Equilibrium 
(64] 

19 |Rezaice-Pajand [Rectangular [5 10 Complete | Complete | Compatibility[no [Can be used as 
and Yaghoobi linear linear Equilibrium transitional ele- 
[64] ment 

20 | Rebiai and Be- | Quadrilateral [4 12 [Incomplete [Complete [ Compatibility! yes |= 
lounar [49] fourth-order _| linear 

2 | Rebiat and Be- [Quadrilateral [4 12 [Incomplete [Complete [ Compatibility! yes |= 
lounar [50] fourth-order _| linear 

2 | Rebiai, Saidani, | Quadrilateral | 1 12 [Incomplete [Complete | Compatibility] yes _|- 
and Bahloul [51 fourth-order | linear 

2B_| Rezaice-Pajand | Quadrilateral [5 W Complete [Complete | Compatibility| yes |= 
and Yaghoobi second-order | second-order | Equilibrium 

w : Quadrilateral [9 18 Complete [Complete [Compatibility] yes |= 
and Yaghoobi second-order | second-order 
(67] 

25 |Hamadi, Ayoub, [Rectangular [5 10 Complete [Complete | Compatibility] no 
and Maalem [31] linear linear Equilibrium 


authors took an advantage of linear strain components and evaluated the 


performance of this element in analysis of both thin and thick planes. 
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One of the first applications of the strain formulation in finite element 
analysis dates back to 1972 when Ashwell and Sabir [7] proposed a rectangu- 
lar cylindrical shell element. Later, in 2004, Djoudi and Bahai [24] developed 
a strain-based shell element, which included openings and cut-outs. They 
utilized this element to perform vibration analysis of shell structures, and 
concluded that the strain-based element is more economic than the conven- 
tional displacement-based elements. In a comparative study, Hamadi et al. 
[32] investigated the performance of strain-based shell elements in comparison 
with the displacement-based elements. This study can be considered an ex- 
tension of their previous research about plate elements [30]. In 2015, Mousa 
and Djoudi [45] performed vibration analysis on circular cylindrical shells 
with oblique ends, by using a new strain-based triangular shell element. For 
this element, they assumed linear curvature and third-order normal in-plane 
strains. The most-recent strain-based shell element is a flat triangular hybrid 
element proposed by Rezaiee-Pajand and Yaghoobi [69]. Trefftz functional 
was used in this element to formulate independent internal and boundary 
fields. It must be noted that this work is not the only available hybrid 
strain-based element and there are other similar studies in the literature. 
For instance, To and Liu [77] also proposed a triangular shell element, based 
on the Hellinger-Reissner hybrid strain formulation. Moreover, the hybrid 
formulation was also utilized for other types of structures. For example, in 
2017, Rezaiee-Pajand and Yaghoobi [68] developed a hybrid plane element 
with assumed strain field. 

In addition to the plate and shell element, the assumed strain approach 
was also utilized for developing three-dimensional finite elements. Belounar 
and Guerraiche [16] formulated a 3D eight-node brick element by assump- 
tion of linear strains. In this formulation, only compatibility criterion was 
imposed to the assumed strain field. Guerraiche, Belounar, and Bouzidi [29] 
proposed another variant of this element by using a different assumed strain 
field. In this new element, they included Poisson’s ratio of the material in 
the strain field of the element. In the most-recent study, Messai, Belounar, 
and Merzouki [44] suggested a nine-node brick element with linear assumed 
strain field. In one of the most recent studies in this field, Rezaiee-Pajand, 
Gharaei-Moghaddam, and Ramezani [60] utilized strain-based formulation 
and developed a cracked plain element for analysis of cracked structures with 
open stable cracks. 


7 Discussion and conclusion 


In this article, which is the first part of a two-part study, the basic formula- 
tion steps of the assumed strain approach for developing plane elements were 
presented. In addition, twenty-five of the available strain-based membrane 
elements were reviewed, using the unified notations. The reviewed elements 
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were categorized into two groups of triangular and quadrilateral elements. 
According to the available literature, most of the accessible elements were 
formulated by using strain fields, which do not satisfy the equilibrium equa- 
tions. Moreover, in many of these elements, the assumed strain field seems 
to be resulted from the trial-and-error process. Because no clear justification 
was provided by the authors for their selections. Despite this fact, in some 
cases, especially the elements proposed by Rezaiee-Pajand and his colleagues, 
the concept of Taylor expansion was considered in the selection of assumed 
strain components, and in many of these elements, both compatibility and 
equilibrium criteria were imposed. It is known for the researchers that utiliz- 
ing incomplete polynomials for strain field might lead to incapability of the 
element to include Poisson’s effect of strain states. Therefore, in some cases 
the authors tried to consider those strain fields with complete polynomial 
approximation. It is worth mentioning that the other elements, which do not 
consider this limitation, can also provide accurate results. 

According to the presented review and in some cases, the same elements 
were proposed by different authors, independently. Therefore, this is one of 
the main motivations of the present study to provide a comprehensive ref- 
erence to be used by authors, which helps prevent such duplicate element 
formulation. In addition, it seems that most of the researchers who worked 
in this field only followed their own research line and did not take advantage 
of the other’s experiences. Accordingly, there are many unanswered ques- 
tions about the performance of the strain-based elements that remain to be 
answered. For instance, it is widely known that some patterns for distri- 
bution of the element degrees of freedom resulted in the singular geometric 
matrix, which was shown by A. In addition, it was concluded that distribu- 
tion of the degrees of freedom between the element nodes had a considerable 
effect on the element performance. Despite these indications, it is not clear 
how to prevent this problem and choose an optimal configuration. Due to 
these warnings, the application of the higher-order strain fields for membrane 
elements requires further investigation. According to these points, in the sec- 
ond part of this study, extensive numerical investigations will compare the 
performance of the reviewed membrane elements. 
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